Critical Temperature and Condensate Fraction of a Fermion Pair Condensate 
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We report on measurements of the critical temperature and the temperature dependence of the 
condensate fraction for a fermion pair condensate of 6 Li atoms. The Bragg spectroscopy is employed 
to determine the critical temperature and the condensate fraction after a fast magnetic field ramp to 
the molecular side of the Feshbach resonance. Our measurements reveal the level-off of the critical 
temperature and the limiting behavior of condensate fraction near the unitarity limit. 



Ultracold fermionic atoms endowed with tunable in- 
teraction offer an ideal testing ground for many-body 
theory [l], 0, 0|. Near a Feshbach resonance, the system 
provides an access to Bose- Einstein condensation (BEC) 
and Bardeen-Cooper-Schrieffer (BCS) crossover [HH, 0,0, 
H, Q . Investigation of the superfluid transition tempera- 
ture T c in the BCS-BEC crossover regime is expected to 
unravel the underlying physics in a strongly interacting 

system^ M M M, M El M E3 ■ 

Measuring T c in the strongly interacting regime poses 
two formidable challenges: thermometry and identifica- 
tion of the emergence of a molecular condensate. Unlike 
the weakly interacting regime, the in-trap cloud size is 
not directly related to temperature because the cloud 
is strongly distorted by the interaction. The onset of 
the BEC transition near the unitary regime is also dif- 
ficult to identify because the bimodality of the distribu- 
tion in real space is smeared out by the strong interac- 
tions. The time-of-flight (TOF) technique is not applica- 
ble either because the interaction e nerg y is converted to 
kinetic energy during expansion @, llcf . A widely used 
method of temperature measurement is to sweep the 
magnetic field isentropically to a nonintcracting value 
at which the temperature is deduced from the density 
profiled, [H, US El- However, this method is not ap- 
plicable to the BEC side of the resonance because of the 
short molecular lifetime [22l. |23|. 

To observe the momentum distribution of fermion 
pairs, JILA and MIT groups 0, i] utilized a rapid 
magnetic-field ramp to convert fermionic pairs to tightly- 
bound molecules while preserving their center-of-mass 
(COM) momentum. This process is referred to as pro- 
jection. Since molecules after the projection are weakly 
interacting, it is possible to measure the condensate frac- 
tion and the COM momentum distribution of pairs by 
using a TOF technique. 

In this Letter, we report on the measurements of T c 
and the condensate fraction of fermion pairs of 6 Li in the 
BCS-BEC crossover regime using the magnetic field ramp 
technique. We determined the temperature from the 
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COM momentum distribution of fermion pairs by carving 
out a slice of the momentum distribution using the Bragg 
diffraction [H, [2^| . Our thermometry based on the Bragg 
diffraction does not suffer from the distortion of a molec- 
ular cloud due to inhomogeneity of the magnetic field 
pulse for molecular dissociation, since the Bragg pulse 
is applied before the dissociation field pulse. With the 
Bragg spectroscopy, we are also able to separate the zero 
momentum component from the rest of the cloud, and 
unambiguously identified the emergence of a molecular 
condensate. In fact, the critical point of BEC can be 
identified with the point at which the number of Br agg - 
diffracted molecules begins to increase precipitously [26j. 
The measured dependence of T c on the s-wave scattering 
length can be explained by theory of weakly interact- 
ing bosons in the weakly interacting regime, whereas in 
the strongly interacting regime, T c deviates significantly 
from the theoretical prediction and levels off near the 
unitarity limit. Concurrently, the measured dependence 
of the condensate fraction on temperature shows limiting 
behavior. 

In our experiment, we employed an all-optical creation 
of Li pairs from atoms in the hyperfine ground states of 
\F,m F ) = |l/2,l/2) (= |1» and \F,m F ) = |l/2, -1/2) 
(= 1 2)). We captured atoms in a cavity-enhanced op- 
tical dipole trap, which were loaded directly from a 
magneto-optical trap27j. A cavity-enhanced 1064 nm 
laser achieved a trap depth of /cb x 2 mK with a beam 
waist of 260 /.im. We then transferred the atoms into a 
focused single-beam optical trap with a waist of 27 /mi. 
The radio-frequency field was applied to produce equal 
populations in the |1) and |2) states. The evaporative 
cooling was initially performed at 834 G, and then the 
field was adiabatically ramped to the magnetic field J5 evap 
where fermion pairs were produced by the final evapora- 
tion. 

The temperature was controlled by tuning the final 
trap depth of the optical trap in the evaporation. The 
number of created molecules ranges from 1.0 x 10 5 to 
1.0 x 10 6 , depending on the evaporation field and final 
trap depth. Trap frequencies are w ra( j/27r = 90.1\/PHz 
and u>ax/27r = V0.57 P + 0.33 B Hz in the radial and ax- 
ial directions, respectively, where P is the laser power of 
the optical trap in mW and B is the strength of mag- 
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netic field in G. The laser power ranges from 30 mW to 
170 mW, which corresponds to the trap aspect ratio from 
30 to 60 under the experimental condition. 

The time sequence of the Bragg diffraction spec- 
troscopy is shown in Fig. [1] After evaporative cooling 
at various magnetic field -B cva p, we held atoms/molecules 
for 100 ms to damp out possible excitations and then we 
turned off the optical trap (t = 0). Simultaneously, we 
turned off the magnetic field rapidly at a sweep rate of 
15 G///S so that the interaction will not be converted to 
kinetic energy during expansion, and therefore the follow- 
ing Bragg spectroscopy does not suffer from resonance 
shift and the broadening (24j. Since the ramp time of 
magnetic field is much shorter than the collisional time of 
molecules, the growth time of the condensate, and other 
time scales of the dynamics, the velocity distribution of 
molecules after the ramp should reflect the initial COM 
motion of fermion pairs. After switch-off of the optical 
trap, the Bragg pulse is applied to the falling molecular 
cloud along the axial direction of the trap(28[. The num- 
ber of molecules, which are momentum-selected by the 
Bragg diffraction, is counted for each frequency difference 
between the two Bragg beams. In order to detect the 
molecules, we ramped up the magnetic field across the 
Feshbach resonance (834 G) to dissociate the molecules. 
Then, we again switched off the magnetic field to cross 
the Feshbach resonance nonadiabatically before taking 
images. At this stage, the atomic cloud has already been 
expanded sufficiently, and therefore the re-association of 
molecules is negligible. 

Figures [2] (a,b) show the density profile of molecules 
created at 780 G at slightly below T c without (a) and with 
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FIG. 1: Time sequence of the experiment. After evaporative 
cooling at varying magnetic field B Bva . p , the optical trap is 
tuned off at t = 0. At the same time, the magnetic field is 
rapidly turned off to convert fermion pairs to tightly-bound 
molecules. The Bragg pulse is applied at t = 500 fjs. Af- 
ter a 3-ms free fall, we ramped up the magnetic field across 
the Feshbach resonance to dissociate the molecules. Then, 
the magnetic field is again switched off to cross the Feshbach 
resonance nonadiabatically before taking absorption images. 
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FIG. 2: Bragg diffraction spectroscopy. Molecules slightly 
below T c without (a) and with (b) the zero momentum com- 
ponent that is diffracted by the Bragg pulse. The profiles 
represent the radially integrated density profiles and the in- 
sets show the two-dimensional images, (c) The number of 
diffracted molecules at the Bragg resonance. T c is identified 
with the point (indicated by the arrow) at which the slope 
changes, (d) The Bragg spectrum of strongly interacting 
molecules at 780 G. The condensate fraction and tempera- 
ture are deduced from this spectrum (see text). 



(b) application of the Bragg pulse at the resonance for 
zero momentum molecules. The Bragg diffraction condi- 
tion is described as hAv = (2hk) 2 /2m m + Hq- (2ftk)/m m ; 
here, Av is the frequency difference between the Bragg 
beams, m m is the mass of a molecule, k is the wave 
vector of the Bragg beam, and q is the initial wave 
vector of diffracted molecules. For molecules at rest, 
hAv = (2fifc) 2 /2m m = h x 147.75 kHz is the Bragg reso- 
nant condition. Although the sharp top of the profile in 
Fig. [Ha) hints the existence of the condensate, it is not 
straightforward to pinpoint the emergence of the conden- 
sate. Using the Bragg diffraction, a small condensate can 
be clearly detected as shown in Fig. EJb). By counting 
the number of diffracted molecules while changing the fi- 
nal trap depth, we have observed a sudden increase in the 
number of diffracted molecules (Fig. [He)). This sudden 
increase of diffracted molecules shows the onset of molec- 
ular condensate as indicated by the arrow in Fig. HJc). 
Once T c is identified, we fix the final trap depth and mea- 
sure the number of diffracted molecules by varying Av 
to obtain the momentum distribution. 

Figure [D(d) shows a typical Bragg spectrum below 
T c . The narrow peak at the center shows the condensed 
molecules, and the peak width arises solely from the time 
duration of the Bragg pulse. No resonance shift or broad- 
ening due to interaction was observed in our experiment, 
indicating that the inter-molecular interaction is negligi- 
ble at zero magnetic field. The condensate fraction and 
the temperature are deduced by fitting the data with a 
function which is obtained by the convolution of a bi- 
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FIG. 3: Temperature dependence of the condensate fraction. 
The condensate fraction is plotted against T/Tf, where Tp 
is the Fermi temperature of a harmonically trapped system. 
Dotted lines show guides to the eye. The solid curve shows the 
condensate fraction of noninteracting bosons in a harmonic 
trap. The vertical error bars involve statistical and systematic 
errors, the latter arising from possible inelastic collisions dur- 
ing TOF after proiectionpj^]. The inset shows the condensate 
fraction at 760G obtained from the Bragg spectroscopy (filled 
squares) and from fitting the profile after projection with a 
Thomas-Fermi distribution for the condensate and a Gaus- 
sian distribution for the thermal component (open squares). 
The two sets of data agree within experimental uncertainty. 



modal distribution and Rabi oscillations, 



molecules observed after projection, hp is the Fermi wave 
number, and a is the atomic s-wave scattering length [2^|. 
The temperature is plotted in units of the Fermi temper- 
ature Tp = fowhot^-ZV) 1 / 3 /fee of a harmonically trapped 
noninteracting system with Who being the geometrical av- 
erage of the trap frequencies. The number of atoms N 
is measured from absorption images at high magnetic 
field [30J]. Dotted lines in Fig. [3] show guides to the eye 
produced by fitting the five lowest data points with a 
linear function. Our data show that as we approach the 
Feshbach resonance, the condensate fraction curve shifts 
to the low temperature side and eventually approach the 
limiting curve (see data points for 760G and 780G). To 
compare our Bragg spectroscopy and the conventional bi- 
modal fitting scheme, we plot in the inset the condensate 
fraction for 760 G determined with the Bragg diffrac- 
tion method (filled squares) and the condensate fraction 
determined using the bimodal fitting to the profile after 
the projection and the ballistic expansion (open squares). 
For the bimodal fitting, we use a Gaussian function for 
the thermal component and a Thomas-Fermi profile for 
the condensate. We can clearly see that our Bragg spec- 
troscopy and a bimodal fitting give the same condensate 
fraction. 

Next, we focus on the critical temperature which can 
be used as the acid test of crossover theory. Figure H] 
shows the measured critical temperature T c plotted in 
units of Tp using the BE distribution (solid circles) , and 
T c determined using the Maxwell-Boltzmann (MB) dis- 
tribution (open circles). Since it is difficult to judge 
whether the BE distribution fits the data better than 
the MB distribution, we also use MB distribution by 
replacing the g 5 / 2 term of Eq. [T] with its argument 
to determine the temperature. Assuming the molecu- 
lar scattering length is 0.6 times the atomic scattering 
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where A, B and T are fitting parameters, g 5 / 2 (z) is a 
polylogarithm function defined by g n (z) = YlkLi z k /k n , 
and Av q — Av — 147.75 kHz. v is equal to (2r 7T )^ 1 , 
where r n is the time duration of the Bragg pulse for the 
7r-pulse condition. The term with the delta function cor- 
responds to a condensate, and the g 5 / 2 term to a thermal 
component described by the Bose-Einstein (BE) distribu- 
tion. Here we use the momentum distribution function 
of noninteracting bosons to describe the thermal compo- 
nent, and therefore the chemical potential is taken to be 
zero. 

Figure [3] shows the temperature dependence of molec- 
ular condensate fraction, N /N m , measured at 690, 720, 
740 G, 760 and 780 G, which correspond to (fcpa)" 1 ~ 
3.4, 1.6, 0.9, 0.6, 0.4, respectively. Here N is the num- 
ber of condensed molecules, N m is the total number of 
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FIG. 4: T c /Tf versus dimensionless interaction parameter 
(fcFa) -1 . Closed circles show T c determined by the BE dis- 
tribution, and open circles show T c determined by the MB 
distribution. The solid curve shows the theoretical prediction 
of weakly interacting bosons in a harmonic trap[3l|. where we 
assume that the molecular scattering length is given by 0.6 
times the atomic scattering length [22|. 
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length |22l| . the theoretical prediction of the weakly inter- 
acting bosons shows good ag reement with the measured 
T c /Ty for (fepa) -1 > 1.5[3l|. In the strongly interact- 
ing regime, T c /Tp gradually shift downward, and even- 
tually levels off for (fcpa)" 1 < 0.6 at T C /T F = 0.24(0.03) 
from BE distribution and 0.19(0.03) from MB distribu- 
tion. The measured T c /T-p near unitarity is close to 
a theoretical prediction of T c /Tp = 0.29[19( and previ- 
ously reported experimental results of 0.27(0.02) [Tj| and 
O.29(+O.O3/-Q.O2)0, both of which were deduced from 
the Thomas- Fermi profile of the gas. It is remarkable that 
within experimental uncertainties T c /Tp remains con- 
stant for (/cfM)" 1 < 0.6. At the unitarity limit our result 
of Tc/T-p is consistent with the theoretical predictions [l9j 
and the previous measurements (l9l. l2lj. indicating that 
our thermometry is valid at the unitarity limit. Our find- 
ing of the limiting behavior may suggest that the univer- 



sal thermodynamics holds true not just at the unitarity 
limit but over an extended region on both sides of the 
Feshbach resonance. 

In conclusion, we used the Bragg spectroscopy to mea- 
sure the critical temperature and the condensate fraction 
of fermion pair condensates in the BCS-BEC crossover re- 
gion. We have succeeded in extracting a condensate from 
a thermal component in the strongly interacting regime 
with high sensitivity and applied the technique to study 
the thermodynamics of strongly interacting fermions. 
Our findings of the limiting behavior in the condensate 
transition temperature and the condensate fraction curve 
near the unitarity awaits deeper understanding of the 
BCS-BEC crossover physics. 

The authors acknowledge S. Inouye and M. Kozuma 
for comments and discussions. 
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